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1. Introduction 

Stars (including the Sun) exert not only gravitation, but also radiation pressure on 
nearby bodies. At the same time, it is well known that dust particles are characterized by a 
considerable sailing capacity (cross-section to mass ratio), and hence are subject to a sizable 
effect of light pressure from the star, being one of the possible mechanisms for the formation 
and evolution of gas-dust clouds. 

In the classical planar restricted circular three-body problem, two large masses, mi and 
77i2, rotate in planar circular Keplerian orbits, while a third particle of negligible mass 
moves in the same plane as the two larger bodies under their gravitational pull. However 
the classical model of the restricted three-body problem is not valid for studying the motion 
of material points in the solar system where the third mass has considerable sailing capacity 
(for example cosmic dust, stellar wind, etc). Thus it is reasonable to modify the classical 
model by superposing a radiative repulsion field, whose source coincides with the source of 
the gravitational field (the Sun), onto the gravitational field of the main body. 

This problem was called the photo- gravitational restricted three-body problem by Radzicv- 
skii [ill . In a later work [2(| , Radzievskii performed a complete treatment of the behavior of 
the equilibrium points. In both papers, however, Radzievskii, who was primarily interested 
in the solar problem, only treated a limited range of radiation pressure parameters (in 
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particular when only one massive body is luminous) and did not consider the question of 
the linear stability of the equilibrium points. 

The photogravitational restricted three-body problem has been treated by several au- 
thors. In 1970 Chernikov [H[ investigated the stability of the collinear equilibrium points L\, 
L2, and L3, as well as the L4 and L5 points, and discussed the modifications of the results 
brought about through the Poynting-Robertson effect, but again only for the Sun-planet 
problem. Later, in 1985, Simmons et al. [23[ investigated the existence and linear stability 
of the libration points, while, in 1996, Khasan Q studied librational solutions of the pho- 
togravitational restricted three body problem by considering both primaries as radiating 
bodies. 

Nonlinear stability of the triangular libration points was investigated by Kumar and 
Choudhry [10], who extended the work of Radzievskii by analyzing the stability of the 
triangular points for all values of the parameters which describe the radiating effects of the 
primaries. They found that, except for some cases, the motion is stable for all values of 



the radiation reduction factors and for all values of \i < 0.0285954 . They also studied [11 1 
the stability of L4 and L5 under the resonance conditions u\ = 2ui2 and u>\ — 3^2- Later, 
Gozdziewski et al. 0] also studied the nonlinear stability of the triangular libration points. 
Their study of the stability of the libration points when the resonances do not exist, and in 
the fourth order resonance case, shows that for some values of parameters these points are 
stable, and for others they are unstable. In the case when the third order resonances exist, 
the triangular libration points are always unstable. 

Kunitsyn and Polyakhova [l2j gave a retrospective review of many aspects of the libration 
point positions and their stability for all values of radiation pressure and mass ratios. 

In this paper we investigate the Lyapunov stability of the triangular solutions for all 
possible values of the reduction coefficients given by n\ and K2, which represent the ratios of 
the difference between the gravitational force and the radiative force of the bodies of mass 
mi and ui2 , respectively. 

We would like to stress here that some results presented in this paper have been obtained 
in the above references. The results of this paper should be viewed as a step towards 
describing the dynamics of the photogravitational restricted three-body problem. We note 
that, seen as a whole, the problem displays a number of interesting features that are not 
apparent in previous treatments. 

The paper is organized as follows: section [2] contains preliminaries, where we recall the 
earlier results on the photogravitational planar restricted three body problem, and introduce 
the Hamiltonian function and the equilibria of the system, along with the condition on the 
masses and radiation pressure values for the stability of the linearized problem. 

In Section the study of the existence and linear stability of equilateral libration points 
is presented. 

Finally, in Section |4l we apply the Arnold- Moser theorem to examine the condition of 
non-linear stability, excluding the resonance cases up to the fourth order. We obtain a 
condition for stability in the parameter space of K\, K2 and \x which may be expressed as 
a quintic in /x, plot the relevant surface, and show that this polynomial reproduces the 
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classical result of Deprit and Dcprit-Bartholome [6| when K\ = K2 = 1. Special cases when 
this polynomial reduces to a quartic (and therefore explicit algebraic expressions for the 
surface may be obtained) are determined. 

2. Hamiltonian formulation 

We consider an infinitesimal mass particle moving in the photogravitational field of two 
masses, termed the primaries, mi and m.2, with both masses in circular orbits around their 
common center of mass. The two primaries are sources of radiation, with the parameters 
Ki and K2 characterizing the radiation effect of mi and m-2 respectively. Similar to the 
planar classical case, the motion of the infinitesimal mass takes place in the same plane of 
the primaries. 

We should note that, in contrast to the classical restricted three-body problem, in the 
photogravitational problem the force acting on the particle depends not only on the param- 
eters of the stars (temperature, size, density, etc.) but also on the parameters of the particle 
itself (size, density, etc.) 

The photogravitational version of the restricted problem presented here is derived in a 
similar way to the classical problem (see [26jj): in fact we start from a presentation and 
notation very similar to that used by Simmons (23|. Units are chosen so that the unit of 
mass is equal to the sum of the primary masses, mi + 7712 = 1, the unit of length is equal 
to their separation, and the unit of time is such that the angular velocity oj = 1. We also 
set the gravitational constant G = 1. For definiteness we also take fj, = to 2 /( to i + m ?)-, 
< /1 < 1 so that mi = 1 — fi and m 2 = H- 

We have fixed the center of mass at (0, 0) and the primaries, toi and m,2, at (— /i, 0) and 
(1 — /i, 0), respectively. The forces experienced by a test particle in the coordinate system 
rotating with uj = 1 and origin at the center of mass are then derivable from the potential 

U { x,y) = Ux 2 +y 2 )+ a ^± + ^ (1) 
2 n r 2 

where {x, y) are the coordinates of the test particle, 

r\ = (x + n) 2 +y 2 , rl = (l-x-fi) 2 +y 2 

are the distances from the masses mi and m.2, respectively, and a, ft represent the effects 
of the radiation pressure from the two primaries. 

The Jacobi constant of the problem is given by 

Cj(x t y, x, y) = \(x 2 + y 2 ) + \{x 2 + y 2 ) + + ^. (2) 

2 2 ri r 2 

After introducing the canonical coordinate system (x,y,p x ,p y ) 

dx dy 



3 



one obtains the Hamiltonian function 



H(x,y,p x ,py) = -(p x +p v ) +p x y~p y x + — h — . (3) 

In this canonical formulation the problem of determining the libration points consists 
of finding all real solutions of the system of four algebraic equations given by setting the 
Hamiltonian equations of motion equal to zero, that is 

Px + y = 0, p y - x = 0, 

a{l - fi)(x + fj.) Pn(x-l + n) 

-Py + zs + ~3 = °> ( 4 ) 

' 1 '2 

For later convenience we introduce the parameters nf = a and k?, = /3 as used by 



Schuerman [221] ■ Defining 

we find that the solution of Q for the coordinates of the triangular libration points L4 and 
L 5 is (see [10j|): 



(6) 



t/L 4 = -yi B = KiK 2 Vb, p VLi = p VL5 = X L4 . 

This implies that 

ri = a x / 3 = Ki, ra^^^^, (7) 

showing that a necessary condition for the existence of triangular points is k% > and 
K2 > (see In fact, as shown by Schuerman the points lie at the intersections of 

the circles defined by 0, and so they exist provided further that K\ + K2 > 1- We remark 
that the existence of the equilibrium solutions is critically governed by the numerical value 
of b. 

These points are the vertices of two triangles, of sides k%, K2 and 1, based on the line 
joining the primaries, and so L4 and L5 are known as triangular (Lagrange) libration points. 
If we set a = 1 and fa = 1, the restricted photogravitational three body problem is reduced 
to the classical case. We consider here only the stability of L4, however all conclusions about 
the stability of L4 can be extended to L$ just applying the symmetries of the photogravita- 
tional problem, namely 

{x,y,p x ,p y ,t) i-> (x,-y,-p x ,p y ,t). (8) 



4 



3. Linear stability 



To investigate the stability we use Birkhoff's procedure |3| for normalizing the Hamilto- 
nian in a neighborhood of the libration point. We start the study of stability by finding the 
first order variational equations, which we then use to determine the eigenvalues. By the 



Arnold- Moser theorem (see [16]), it is known that a necessary condition for stability of the 



L4 point is that all eigenvalues should be pure imaginary. 

Using a linear canonical transformation, we first shift the origin of the coordinate system 
to the L4 point and expand the Hamiltonian in a power series of the coordinates. To this 
end we define new coordinates q\ , 92 , Pi, Vi by 

x = x Li +q 1 , Px=Px Li +Pi, 

y = VL±+ 92, Py = Py Li + Pi, 

in terms of which the equations of motion are 

dq^_dH_ dp^ _dH_ • — 10 

dt dpi ' dt dqt ' 

In these new variables the solution ^ corresponds to the equilibrium state qj = pj = 0, 
for j = 1, 2.F We now expand the Hamiltonian in a power series of the coordinates, up to 
fourth order. Since the expansion is made in the neighborhood of an equilibrium point, the 
constant term (the value of the Hamiltonian at equilibrium) can be neglected, and the linear 
part must vanish. The expanded Hamiltonian can be written as 

00 

H = Y J H 1 , (9) 

3=0 

where Hj are homogeneous polynomials of degree j in the new variables. Calculating this 
expansion to fourth order we find 

H = H{x Li ,y Li ,p XLi ,p yL J, 
Hi = 0, 

H 2 = \{p\+pl) + q2Pi - QiP2 + {A-\)q\ + B qi q 2 ~{A+\) q\, 

Hz = /l3000 Qi + ^2100 Q1I2 + ^1200 92 91 + ^0300 ?2i 

Ha = /14000 Qi + ^3100 9i92 + ^-2200 qlli + ^1300 Q1Q2 + ^0400 Q2, 

where 

A = -\ + y [>«? + (1 - A*)*!] , 

B EE -|^-[(l- M ) K 2 (K 2 + 1 _ /t 2 )+MK 2 (K 2 + 1 _ K 2 ) ]_ (10) 

The values of A and B used in Hi are conveniently defined for simplifying expressions which 
appear later. The coefficients of H3 and H4, calculated using Maple, coincide with those 
given in Kumar and Choudhry [10( and are listed in Appendix A. 
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The system of linear differential equations derived from the quadratic term, H 2 , describes 
the tangent flow around L4. Linear stability is determined by the character of the associated 
eigenvalues, which are the roots of the characteristic equation 

A 4 + A 2 + 9^(1 — fi)b= 0, (ff) 

whose eigenvalues are found to be 



Ai 2 = ±1 



^3,4 



1 


v/1 


- 36/i(l 


- M )6 


2 H 




2 




1 




- 36/x(l 


-fx)b 


2 " 




2 





(12) 



For linear stability all of these eigenvalues should be pure imaginary, which is the case if 

< 36/i(l - fi)b < 1. (13) 

Substituting (O in (fT3")) , we obtain 

k|-1\ 2 1 + 4 - 1\ 2 ^ -36/i 2 + 36/i- 1 



1 - x ' a — - < — f — - => 1 " — - > f ; \ . (14) 

\ 2kiK2 / 36/i(l — /i) \ 2ki«2 / 36/i(l — /i) 

Since the left-hand side can be zero, we require 

-36u 2 + 36u - 1 9 

' < ^ 36ju 2 - 36// + 1 > 0, (15) 



36/i(l — /i) 



as /i € [0, 1]. This equation is satisfied for /1 < | — and /1 > i + and so, in these 
subintervals, the eigenvalues (|12p are distinct and pure imaginary. 

Throughout this paper, we restrict ourselves to the case 

1 /2 

H < - - \- 0.0285954 = /i* 

for which the linear stability conditions are fulfilled for all values of n\ and k 2 - 

Note that Lagrange's classical criterion for stability is recovered in the limit of no radi- 
ation pressure (k\ — > 1, n 2 —> 1, b —> 3/4). In this case the inequality (fT4|) for the existence 
of two pairs of pure imaginary values becomes 

27^(1 — fi) < 1, or equivalently /i(l — fi) < 1/27. 

We remark that this result implies the stability of the L4 point in the classical planar 
circular restricted three-body problem when the mass ratio parameter satisfies < fi < 
fiR = 5(1 — -^ip) ~ 0.0385208, where ijlr is known as the Routh value. 
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We write the four eigenvalues as Hwi, ±iu} 2 where the strictly positive numbers u>\ and 
u>2 {frequencies) are determined by the set of relations 

u?- A 2 - 1 + M „ 2 - A 2 _!- M ( W ) 



where M = y/T— 36/i(l — Jijb. The expressions (fl6|) show that 

< w 2 < -r= < wi < 1. (17) 
v2 

We remark that u>i and W2 and the coefficients of all Hj are functions of the parameters /j,, 
K\ and K2- 



4. Normal form and non-linear stability 

Having established linear stability, the next step is to transform the Hamiltonian into 
its Birkhoff normal form. This normalization allows us to apply Arnold's theorem [l[ to 
investigate the stability of the L4 point for mass ratios fj, < /x*, except for the resonant 
cases. From the coefficients of the normal form, Arnold's thorem constructs a determinant, 
-D, defined later, which, when non-zero, establishes the stability of the equilibrium point. 
To obtain the Birkhoff normal form we use the Lie series method, with the calculations 
being performed using Maple. To determine the domain of applicability of the results, it is 
necessary to obtain the resonances, which we now calculate. 

4-1- Existence of resonances 

In this section we study the values for which the frequencies at the equilibrium are in 
resonance. Since we are working to fourth order in the normal form, we thus need to consider 
resonances up to fourth order of the triangular libration points. 

The stability analysis of the L4 point can be carried out if the frequencies u)\ , ui2 satisfy 
the non-resonance condition 

C1W1 + c 2 ui2 7^ (18) 

for all integers c\, C2 such that \c\ \ + |ca| < 4. This condition is violated for ui\ = 2ui2 and 
uji = 3w2, and, obviously when uj\ = u>2- 

We start the discussion by noting that the first-order resonance appears when one of the 
frequencies is zero. In our case, this is possible if either ji = or K\ + K% = 1. For fj, = 
the restricted photogravitational three body problem is reduced to Kepler's problem in a 
rotating coordinate frame. The L4 point in this case is evidently unstable in the sense of 
Lyapunov. 

The second case has no classical equivalent: L4 and L5 coincide with the inner collinear 
point L\ (see Simmons et. al. 23]). In particular if k% = 1 and k\ — > the libration 
points L4 and L5 move from their classical equilateral positions onto the luminous mass and 
coalesce there with the inner libration point L\. It is clear that the motion will be unstable. 
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We have seen above that if Ki = K2 = 1, then fj, = fiR and the triangular libration points 
are stable in the sense of Lyapunov. This agrees with the result of Meyer Placian and 
Yaguas ljj, who demonstrated the stability of the Lagrange equilateral triangle points, L4 
and L5, in the plane circular restricted three-body problem when the mass ratio parameter 
is equal to = [ir, the critical value of Routh. 

We turn to the two remaining cases. The second order resonance wi = 0J2 occurs when 



/i = i ± i — -gij. Since /1 > 0, the resonance appears for parameters K\ and K2 satisfying 
the condition b > 1/9. In the first graph of Figure I4TT1 we have plotted the relevant part (i.e. 
with the negative root) of n(b), together with the line fi = fj,*. We see that (other than at 
the exceptional value 6=1) the resonant values only exist for /i > fi* . 
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Figure 1: Graphs of /i(b) for (from left to right) the 1:1, 1:2 and 1:3 resonances. The horizontal line is 

We use the same criterion adopted in the previous case for determining the masses for 
which the resonance uj\ — 2ui2 occurs. In this way, one obtains fi = ~ — \J \~ > so * na * 
the resonance appears for parameters K\ and «2 satisfying the condition b > 16/225 « 0.071. 
The function /j,(b) for this resonance is shown in the second graph of Figure 14.11 We see 
that n(b) = fi* when b = 16/25 = 0.64. Therefore, in contrast to the results of Kumar and 
Choudhry [10(, who claimed that the resonance u>i — 2u>2 occurs for values from b = 0.65 to 
b = 0.95, we find that it occurs for b € [0.64, 1]. 

Finally, the resonance u>i = 3^2 occurs for fi = | — \J j ~ ■ Since fi > 0, the 
resonance appears for values of K\ and «2 satisfying b > 1/50. In this case fi(b) = \x* for 
b = 9/25 = 0.36. Kumar and Choudhry [l(| stated that this resonance occurs for values 
from b — 0.4 to b = 0.95, but we find that, in fact, it appears for b g [0.36, 1]. 

In Table Q] we show the values of /i for k\ — 1 and several values of K2 for the resonances 
cji = 0J2, = 2o;2 and ui\ — Zuj2- Note that when one of the radiation coefficients is 1, and 
there is a 1:1 resonance, the points L4 and L5 are linearly stable with the condition that 
mass parameter /1 is at most 0.0385209. Our result refines that of Kumar and Choudhry [l(| 
and has an influence on the analysis of the non-linear stability of L4, and consequently on 
that of L5. 
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K 2 




oji = 2a>2 


oji = 3a>2 


1 


0.0385209 


0.0349233 


0.0321444 


0.9 


0.0361369 


0.0327702 


0.0301682 


0.8 


0.0342412 


0.0310571 


0.0285955 


0.7 


0.0327266 


0.029688 


0.0273381 


0.6 


0.0315184 


0.0285955 


0.0263345 


0.5 


0.0305638 


0.027732 


0.0255412 


0.4 


0.0298247 


0.0270634 


0.0249269 


0.3 


0.0292741 


0.0265653 


0.0244692 


0.2 


0.0288932 


0.0262206 


0.0241524 


0.1 


0.0286693 


0.0260181 


0.0239662 



Table 1: Values of fj, for k\ = 1. 

4-2. Birkhoff's Normal Form 

In section [3] we determined the conditions under which the eigenvalues are purely imagi- 
nary, thus guaranteeing that the phogravitational Hamiltonian @ satisfies the first condition 
of the Arnold-Moser theorem. The next step is the normalization procedure that transforms 
Hamiltonian (J9j) to Birkhoff's normal form up to fourth order, excluding regions close to the 
resonances established in the previous section. 

To this end, first we make a further change of variables, 
to the Hamiltonian ([3]) that diagonalizes H%: 

m 
3=3 

We use a negative sign on the second term to emphasize the non-positive character of the 
photogravitational Hamiltonian. 

We note that H2 is already in real normal form since it depends only on the actions 
= ^ + fj%, v = 1,2. Also, recall that oji and 0J2 are the linear frequencies analyzed in 
the previous sections. The terms of order greater than two depend on (£1, ^> 771 , 772), but 
not necessarily in such a simple way as 82- The purpose of the Birkhoff normal form is to 
find a near identity symplectic coordinate change 

i-\-j-\-k-\-l— m 

z< m )=z+ £ a ijkl mifj l [fj l 2 +O m+ i (20) 

i+j+k+l=2 

so that the even order terms in the Hamiltonian depend on the new variables z = (£1, , 772) 
in such a way that they can be grouped as powers of Ri and i?2- 

It is implicit in this last statement that the Hamiltonian system's natural resonances 
have been taken into account. Let S be some system of resonance relations, i.e. relations of 



9 



the form X s — Y^i=i m i^i with rrii non-negative integers not all zero. From the viewpoint of 
the general theory of dynamical systems, equilibria and periodic trajectories in Hamiltonian 
systems are all resonant. Indeed, for equilibria, if Ai is an eigenvalue of an equilibrium 
of a Hamiltonian system, then A2 = — Ai is also an eigenvalue. So an infinite number of 
resonance relations of the form 

A s = X s + fc(Ai + A 2 ), k e Z 

are satisfied. This means that we can find a non-divergent generating function (and con- 
sequent canonical transformation) that transforms away all terms of odd degree of the 
Hamiltonian ©. We have, necessarily, to retain part of the even terms corresponding to 
the resonant monomials that appear in This will become clearer in the brief presen- 
tation of the algorithm we have used to implement the Birkhoff normal form. Of course, 
a truncated Birkhoff normal form is a polynomial Hamiltonian that is formally integrable, 
since it is expressed only in terms of the actions Rj, j = 1, . . .,n. For more theoretical 
details we refer the reader to Q and [13] , and for practical implementations, the works 
of [3] and 25 1 . In this paragraph we have tried to stress a point that is usually not made 



clearly when calculating the Birkhoff normal forms for the photogravitational problem: the 
difference between the natural resonances of the 2n complex eigenvalues, and the resonances 
of the real eigenfrequencies uj\ and u> 2 . 

Perhaps due to Hamilton- Jacobi tradition, some engineers and physicists, for example 
Kumar and Choudhry [10j , prefer to work with real variables and mixed generating functions 
which depend on both old and new variables. In this work we use the Lie derivative approach, 
which allows us to use a generating function dependent only on new the variables, with 
the advantage of turning the cumbersome passage to the normal formal more transparent. 
Moreover, with this approach, when the coefficients are not algebraic, but numerical, one 



can write fast codes to construct very high order normal forms (see for example 24|, [25| ) 



Before describing the algorithm and the form of the generating function, a further step 
is necessary to prepare the Hamiltonian for the process of normalization. It is easier to 
perform the manipulations in the complex field and, to this end, we shall make use of the 
change of variables, Rc, and its inverse, Cr: 

R c : Xj = -^{Ij - ifjj), yj = 73 + %)> 
C R : $j = -^(xj + i Vj ) , % = 75 (izj + Vj ) , 
where i = to transform £2! fj 2 ) into the complex Hamiltonian 

H(x 1 , x 2 ,yi, 2/2) = iwiiEiyi +iu 2 X2V2 + ^2 h klk2 i lh x^x^ViV 1 ? ■ (21) 

fel+fe2+il+/2=3 

In fact, the complexification can be done directly from Hamiltonian <|3j) to (f2Tj) . as is the 
case in this work (Appendix B). 

Now, we introduce the generating function G(x\, x 2 , yi, J/2) as a series expansion in 
homogeneous polynomials of degree m > 3 

TCI 

G m (xi,x 2 ,yi,y 2 ) = ^2 G j> where G 7 = 9kik 2 hh x^x^y^y 1 ^, (22) 

j=3 k 1 +k 2 +h+l2=j 
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and the associated canonical transformation Tg, such that 

T G H(xx,X2,y 1 ,y2) = Z(X 1 ,X 2 ,Y 1 ,Y 2 ) = H(T G (x 1 ,x 2 ,y 1 ,y 2 )), 

with 

Z(X 1 ,X 2 ,Y 1 ,Y 2 ) =H + {H,G} + ±{{H, G], G} + ^{{{H, G}, G}, G} H , (23) 

where {H, G} = L^G is the usual Poisson bracket (or Lie derivative) of the functions H and 
G. This means that Z is the time one flow of the flux generated by the canonical system of 
equations associated with the generating function G. To find the Birkhoff normal form and 
the corresponding change of variables (and its inverse) one has to determine the generating 
function that gives the prescribed form discussed above. Collecting powers of (|23jl up to 
order four gives 

Z 2 =H 2 , 

Z 3 =H 3 + {H 2 ,G 3 }, (24) 
Z 4 =Hi + {H 3 , G 3 } + ^{{H 2 , G 3 }, G 3 } + {H 2 , G 4 }. 



Since we are supposing that there are no resonances besides the natural ones (as explained 
above), all terms of order 3 can be removed from (|23j) by setting 

{H 2 , G 3 } = Z 3 — H 3 . 

Solving for Z 3 = gives 

r\ - \" h klk2h i 2 k ka h . 

where 

(k - 1, A) = i(fci - h)ui + i(fc 2 - 1 2 )uj 2 

are the coefficients of the monomials of the Poisson bracket {H 2 ,G 3 }, which are nonzero 
outside the resonance regions. 

For terms of order four, the solution is not so easy since Z 4 ^ 0, and we have to solve 
simultaneously for G4 and Z4, and so on, up to the required order, m, of the normal form. 
In this case one can use the general result [24| for k > 3: 

{H 2l G k } + Z k = F k , F 3 = H 3 , Z 2 = H 2 (26) 

and m m 

F k = — -{G 2 +m, Zk-m} + J]] T ^^2+m : fe-m-2 (27) 

m— l,...,fc— 3 m— 2 

where 

k—m — 2 

H2+m,k—m — 2 = / y T ^ > H2-\-m,k— j— m— 2 } + (28) 

/v ■ 777 A 
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Applying these formulae we can find and G4, which give us the fourth order terms 
that the Arnold-Moser Theorem requires. The algebraic manipulations described so far were 
performed using the software Maple, but any algebraic software could have been used. Note 
that Z± is a complex Birkhoff normal form but, using the transformation Cr given in (|21[) . 
we obtain finally 

Hr 3 = 0, (29) 
Hr 4 = 6 11 (g + vl) 2 + 6u(& + tiM + Vl) + M£a + 
which is the form required by the Arnold-Moser Theorem. 

According to this theorem we have to check whether the determinant 

D(ki,K 2 ) = —{Six Ll>2 - 28l2 LOiU>2 + $22 u\) (30) 

is nonzero. If D(k,i,k 2 ) 7^ for some pair (ki,k 2 ) then, for this pair, the motion is stable 
in the Lyapunov sense provided uj\ 7^ w 2 , wx 7^ 2w 2 or u± 7^ 3w 2 . 

The complex canonical transformation that takes the Taylor series (|3]) to (|2"Tj) is shown in 
Appendix B while the coefficients of the fourth order normal form are given in Appendix C. 
After a sequence of algebraic manipulations we finally obtain the determinant as a rational 
function, whose numerator is a quintic in /1. Explicitly, defining the product of factors 

II = (kx + 1 + k 2 ) (kx + 1 — k 2 ) (ki — 1 + k 2 ) («i — 1 — k 2 ) , 

we can write D as 

1 5 

D(kx,k 2 ) = 

where the denominator, A, is 

X = 2048 n K i 4 K 2 4 (22511^- 1) - 16ki 2 k 2 2 ) {9U(J,(fi - 1) - kx 2 k 2 2 ) 
and the coefficients, Di, of the quintic in the numerator are given by 

(k 2 — K 2 )(5 K\ 2 — 5 — 8 k,\k 2 + 5 k 2 2 )(5 K\ 2 — 5 + 8 K1K2 + 5 k 2 2 ), 



D5 
12960n 3 



Da 
' 2592n 2 



375 kx 2 - 1500 ki 4 + 2250 Kl 6 - 1500 Kl 8 + 375 ki 10 

-250 k 2 2 + 1000 k 2 4 - 1500 k 2 6 + 1000 k 2 8 - 250 k 2 W 
-350 ki 2 k 2 2 + 500 ki 2 k 2 4 - 1450 ki 2 k 2 6 + 925 ki 2 k 2 8 
+ki 4 (150 k 2 2 + 8 k 2 4 - 1650 k 2 6 ) + 1750 kx 6 k 2 2 + 1900 ki 6 k 2 4 
-1300ki 8 k 2 2 , 
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K2 4 



— |- = 3375 ki 2 - 13500 ki 4 + 20250 ki 6 - 13500 ki 8 + 3375 ki 10 
288II 

-1125 k 2 2 + 4500 k 2 4 - 6750 k 2 6 + 4500 k 2 8 - 1125 k 2 10 
-6300 ki 2 k 2 2 + 7495 k 1 2 k 2 4 - 9590 k 1 2 k 2 6 + 5020 k 1 2 k 2 8 
+4205ki 4 k 2 2 + 144ki 4 k 2 4 - 10065 ki 4 k 2 6 
+ 14990 ki 6 k 2 2 + 14565 ki 6 k 2 4 - 11770 ki 8 k 2 2 

~°l = 1125 - 6750ki 2 + 16875ki 4 - 22500ki 6 + 16875ki 8 - 6750ki 10 
288 Ki A H 

+ 1125 ki 12 - 5400 k 2 2 + 11050 k 2 4 - 12700 k 2 6 + 8925 k 2 8 
-3700k 2 10 + 700k 2 12 + 14890ki 2 k 2 2 - 9518ki 2 k 2 4 - 1634ki 2 k 2 
+6012 Ki 2 k 2 8 - 3000k! 2 k 2 10 - 5560 kA 2 2 - 674 Ki 4 k 2 4 
-3712 ki 4 k 2 6 + 7295 ki 4 k 2 8 - 18660 k^k 2 2 - 14298 
-12850 ki 6 k 2 6 + 13440 ki 8 ^ 2 4 + 21440 kA 2 2 - 6710 ki 

~f 1 „ = 3465 - 13860 ki 2 + 20790 ki 4 - 13860 ki 6 + 3465 ki 8 

32 ki 4 k 2 ^ii 

-2790 k 2 2 - 4780 k 2 4 + 4070 k 2 6 + 35 k 2 8 
-2870 ki 6 k 2 2 - 5196 ki 4 k 2 4 + 2950 ki 4 k 2 2 
+2710 ki 2 k 2 2 + 4566 ki 2 k 2 6 + 15912 ki 2 k 2 4 

- ~ D ° . = 5-20ki 2 + 30ki 4 -20ki 6 + 5ki 8 + 10k 2 2 - 36k 2 4 + 22k 2 6 

512 Ki°K 2 4 

-K 2 8 - 10 Kl 2 K 2 2 + 56 Kl 2 K 2 4 + 22 Kl 2 K 2 6 - 10 Kl 4 K 2 2 

-36ki 4 k 2 4 + 10ki 6 k 2 2 

Since the numerator of D is a polynomial of order 5, it is, in general, impossible to 
write an algebraic expression for the values for which /i(Ki,K 2 ) produces D = 0. However 
it is easy to plot the surface D = 0, which is shown in figure 14.21 for the parameter ranges 
Kt <= [0,1], k 2 € [0,1], n € [0,0.3]. 

This surface seems to be equivalent to that shown by Gozdziewski et. al. Q, who 
use slightly different variables. They, however, do not give an explicit expression for the 
determinant (their D4), while our expression allows a better understanding of the structure 
of this function. For example it is possible to see that, for certain combinations of K\ and k 2 , 
the numerator of D reduces to a quartic, and an explicit algebraic expression for (i(ki, h 2 ) 
on the surface D — can be obtained. In the physical region of the parameter space this 
happens when 

K\ + k 2 — 1, K\ = k 2 , or 5(k 2 + k 2 — 1) + 8kik 2 — 0. 

The second of these is particularly interesting, since it includes the purely gravitational case 
Ki = K2 = 1. For these values we have 

9 13041 ^ - 26082 ^ + 14664 [i 2 - 1623 \i + 16 



D = D C 



64 (675 a« 2 - 675 ^+16) (27 ^ 2 -27/i+l) 
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Figure 2: The surface D(fi, m, K%) = 0. 



with the first positive root at 




1 



2898 



1 



1576995 + 966 V199945 w 0.0109136676 



in agreement with the classical result of Deprit and Dcprit-Bartholome. For illustration we 
show the behavior of -D c (/i) in Figure l4~2l 

5. Conclusion 

One of the motivations of this paper was to bring together various results relating to 
the restricted photogravitational three-body problem scattered about the literature, with 
special emphasis on the stability of the L4 and L 5 libration points in the absence of first 
and second-order resonances. As part of this work, we noticed that, though the fourth- 
order Taylor series expansion given by Kumar and Choudhry is correct, one of their results 
apparently does not reproduce the classical case in the appropriate limits. This motivated 
us to use a different approach, namely the Lie triangle method, to calculate the Birkhoff 
normal form for the Hamiltonian. Along the way, it was found that the conditions for the 
existence of resonances given by Kumar and Choudhry had to be slightly modified. 

In this paper we give explicitly the fourth-order normal form, as well as the complex 
canonical transformation used to prepare the Hamiltonian for the Birkhoff normalization. 
We also provide an explicit expression in terms of a rational function, for the fourth-order 
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Figure 3: The classical determinant, D c , for < fi < 0.04 



determinant, D(/j,, k\, K2), arid show that it reduces to the classical, purely gravitational 
case in the limits k± = K2 1- 

Our algebraic analysis seems to corroborate the numerical fourth-order treatment by 
Gozdziewski et. al. [7|. The algebraic treatment shows that the surface D = is, in 
fact, generated by a polynomial that is fifth-order in fi. Knowing the explicit behavior 
of D means that special cases can be studied, and we find three relations in the physical 
parameter space of K\ and K2 for which the numerator reduces to a quartic. In these cases 
explicit expressions for the surface D = can therefore be obtained. One of the cases, 
K i = K 2, includes the classical non-radiational case and we show how the classical result of 
Deprit and Deprit-Bartholome follows as a special case. 
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Appendix A 

The coefficients of the Taylor series expansion for H 3 and H4. hijki is the coefficient of 
Qi ?2 Pi P2 m the Hamiltonian. 



+ («? + 1 - 4) H4 + 1 - K lf - 12 4] - 

8 K 2 

+ [5(k? + 1 - k 2 ) 2 - 4k 2 ] , 

^1200 = || («? - 1 - k 2 )(56k 2 1) + («? + 1 - k|)(56«2 - 1), 

^-0300 = 7r-^-\/& At(5&K 2 - 3) + -^-Vb (1 - ^)(5bnl - 3), 



8k? 



3*tJ-y(^ + l-«2)M + |(^ + l-«2) 4 



8k§ 



1 ^ ^ 

o 4 _ u / 2 1 2^2 2 1 ( 2 _ i _ 2^4 

JK; 2 ~2~^ Kl K 2J K 2 "+" Y6 



5\/& K 2 . o 2\ 

^3100 =—. 3 (1 - A*)(«l + 1 - K 2 ) 

4 1 



. 2 

+ — ; 



4 «r 



t2200 



3(1 - M ) 



4k 2 



5ko& + 



4k? 



2 \ n i 



4) 
+ 1 



3- 

3 - 



7(4 



4? 



4k? 

7(k?-1-k|) 2 



2\2 



4) 



4k| 
35 



+ 



3/i 

4k| 



5k 2 " 

5 (Ji 1 „2\2 35 L„2/„2 



2\2 



5k?6 + (k? - 1 - ni,f - &k?(k? - 1 - k 2 ) 2 - 1 



4k 



4k* 



'1300 



~4k? 



n 2 Vb (1 - At) (k? + 1 - k|)(3 - 74b) 
+ -^-3 KiVb^ («i - 1 - k 2 )( 3 - 7k^6). 



ft-0400 



(1 - /i) (3 - 30k 2 6 + 35^^ 



8k? 



8k 2 



At (3 - 30k?6 + 35k 4 6 2 ) . 
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Appendix B 



The canonical transformation from the variables (qi,Pi) to (Qi,Pi). 

-1 + i j (B- 2ioJx)Qi + (B + 2icu 1 )P 1 (B - 2iu 2 )Q 2 + (B + 2iuj 2 )P 2 



92 



V2 1 ^wi(2o;2 - l)(-3 + 4A- 2u\) ^u 2 (2uj\ - l)(-3 + 4A - 2w|) 



2^ 



-3_+4A_-2ojI 
wi(2o;f - 1) 



1/2 



(Qi+Pi) 



-3 + 4A - 2w, 



w 2 (2o;| - 1) 



1/2 



(Q 2 + P 2 ) 



Pi = 



-1 + i J (2w? + 2iBo;i + AA - 3)Qi + (2w? - 2iBwi + 4A - 3)Pi 



2^2 



(2w| + 2i£w 2 



Vwi(2^ - l)(-3 + 44- 2^) 
4A - 3)Q 2 + (2w| - 2iBu; 2 + 4A - 3)P 2 



VW 2 ^ - l)(-3 + 4A - 2u^) 



-1 + i J [-2iuf + (4A + l)wi - 2B]Qx + [-2iwf + (AA + l)wi + 2B]P 1 
Pl ~ 2y/2 \ ^^i(2uj\ - l)(-3 + 4A- 2uj\) 

[-2k4 + (4^4 + l)u 2 - 2B]Q 2 + [-2*4 + (4A + 1)^ 2 + 2B]P 2 \ 
y/u 2 (2Lj%-l)(-3 + 4A-2bj'i) J ' 

where A and B arc as defined in (fTU)) , with wi and w 2 as in (fTo]) . 



Appendix C 

Applying the transformation of Appendix B to the Taylor series given in Appendix A and 
applying the normal form procedure described in section 4.2, we obtain the following ex- 
pressions for the coefficients of the fourth-order Birkhoff normal after converting to real 
variables variables (xi, x 2 , yi,y 2 ): 



811 



1 / iHQP 1011 HQP 1110 iHQP 2Q01 HQP Q1 



4 V oj 2 2(-2wi + uj 2 ) 

, iHQP 2001 HQP 0120 3iHQP 2010 HQP 1020 

H \ n ^^2020 

iHQP 0021 HQP 21Q0 _ 3iHQP 003 oHQP 3000 

(UJI + UI2) UJl 

, iHQP 0021 HQP 21 .in 
2(-2o>i-w 2 ) 
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c 1 rr np If iHQP 0012 HQP 1200 iHQP 01O2 HQP uw 

022 = —-rHQJr'iiii + - 



Wi + 2 CJ 2 W 2 

iHQP 0201 HQPion iHQP 2ioHQP W c,2 



+ 



+ 



UJ 2 


UJ\ — 2 UJ2 


iHQPoouHQP^oo 


iHQP 02 i HQP W02 


— U!l — 2U! 2 


— UJl +2UJ2 


iHQP 1020 HQPi 10 i 


iHQP 0021 HQP 2 wo 




2 uji + uj 2 


iHQP mil HQP 20W 


iHQP Q i2 HQP20Qi 


UJl 


2 UJl — LU 2 


iHQP m2 nHQP 2ml 


iHQP 002 iHQP 2 ioo 



— 2Wi+Cl»2 — 2 UJi— UJ2 



>12 



1 . 1 ( iHQP m i 2 HQP 12m iHQP 0W2 HQPm 

-- HLjriiii + - — 

4 2 V uji + 2 uj 2 uj 2 



iHQP^oiHQPwn | iHQP 02W HQP wm 

UJ2 UJl — 2 LU 2 

| *HQP W i 2 HQPi 2m iHQP 02W HQP loO2 

—UJi — 2 UJ2 —UJi +2UJ2 

iHQP 1020 HQPnoi iHQP m2 iHQP 2W0 

UJl 2 UJl + U>2 

HIQP0111HQP2010 iHQPoi2oHQP 2 ooi 

UJl 2 UJl — UJ2 

| iHQPoi2 HQP 200 i | iHQP a021 HQP 2W0 
—2uji+uj 2 —2uji—uj2 

In the complex form we have 

H(Q it Pi) = An PiQl + A22 PlQ\ + A12 P1P2Q1Q2 

where 

An = HQ P2020 + — (HQ P W 2oHQP2ow + HQPoo3oHQP 3 ooo) 

UJl 

^iHQPiuqHQPiqh iHQP 2 ooiHQPoi2o iHQP n0 2iHQP 2 ioo 

U! 2 2 UJl — LJ2 2UJI+UJ2 

4 iHQP 00 i2HQPi2oo 4 iHQP mw HQP wm 



A12 = HQP1111 + 



uji + 2 ui2 uji — 2u>2 



| 4iggP 0021 ggP 210 o | AiHQP i2oHQP 2 ooi 
2 uji + 0J2 2 uji — uj 2 

2i 

H (HQPoni HQP2010 + HQPnoi HQP1020) 

UJl 

+ — (HQP 0W 2 HQPmo + HQP0201 HQP W u) 

UJ2 
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and A 22 can be obtained from An by making the changes ui «-» oj 2 and HQPijki — > 
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